home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / zppcon.z / zppcon
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. ZZZZPPPPPPPPCCCCOOOONNNN((((3333FFFF))))                                                          ZZZZPPPPPPPPCCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZPPCON - estimate the reciprocal of the condition number (in the 1-norm)
  10.      of a complex Hermitian positive definite packed matrix using the Cholesky
  11.      factorization A = U**H*U or A = L*L**H computed by ZPPTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE ZPPCON( UPLO, N, AP, ANORM, RCOND, WORK, RWORK, INFO )
  15.  
  16.          CHARACTER      UPLO
  17.  
  18.          INTEGER        INFO, N
  19.  
  20.          DOUBLE         PRECISION ANORM, RCOND
  21.  
  22.          DOUBLE         PRECISION RWORK( * )
  23.  
  24.          COMPLEX*16     AP( * ), WORK( * )
  25.  
  26. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  27.      ZPPCON estimates the reciprocal of the condition number (in the 1-norm)
  28.      of a complex Hermitian positive definite packed matrix using the Cholesky
  29.      factorization A = U**H*U or A = L*L**H computed by ZPPTRF.
  30.  
  31.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  32.      condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))).
  33.  
  34.  
  35. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  36.      UPLO    (input) CHARACTER*1
  37.              = 'U':  Upper triangle of A is stored;
  38.              = 'L':  Lower triangle of A is stored.
  39.  
  40.      N       (input) INTEGER
  41.              The order of the matrix A.  N >= 0.
  42.  
  43.      AP      (input) COMPLEX*16 array, dimension (N*(N+1)/2)
  44.              The triangular factor U or L from the Cholesky factorization A =
  45.              U**H*U or A = L*L**H, packed columnwise in a linear array.  The
  46.              j-th column of U or L is stored in the array AP as follows:  if
  47.              UPLO = 'U', AP(i + (j-1)*j/2) = U(i,j) for 1<=i<=j; if UPLO =
  48.              'L', AP(i + (j-1)*(2n-j)/2) = L(i,j) for j<=i<=n.
  49.  
  50.      ANORM   (input) DOUBLE PRECISION
  51.              The 1-norm (or infinity-norm) of the Hermitian matrix A.
  52.  
  53.      RCOND   (output) DOUBLE PRECISION
  54.              The reciprocal of the condition number of the matrix A, computed
  55.              as RCOND = 1/(ANORM * AINVNM), where AINVNM is an estimate of the
  56.              1-norm of inv(A) computed in this routine.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZPPPPPPPPCCCCOOOONNNN((((3333FFFF))))                                                          ZZZZPPPPPPPPCCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      WORK    (workspace) COMPLEX*16 array, dimension (2*N)
  75.  
  76.      RWORK   (workspace) DOUBLE PRECISION array, dimension (N)
  77.  
  78.      INFO    (output) INTEGER
  79.              = 0:  successful exit
  80.              < 0:  if INFO = -i, the i-th argument had an illegal value
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.